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Abstract. The discovery of deep spectral features in the X-ray spectrum of 
1E1207. 4-5209 focussed the attention of the astronomical community on this radio- 
quiet NS, making it the most intensively observed INS ever. The harvest of X-ray 
photons, collected mainly by XMM-Newton, unveiled a number of unique charac- 
teristics, raising questions on this source very nature. 
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1. Introduction 

Neutron star atmosphere models predicted 
the presence of absorption features de- 
pending on atmospheric composition, but 
high quality spectra, collected both by 
Chandra and by XMM-Newton, did not 
yield evidence for any such feature (see 
Pavlov et al. 2002a and Becker and 
Aschenbach 2002 for recent reviews). INS 
spectra are well fitted by one or more black- 
body curves with, possibly a power law 
contribution at higher energies, but with 
no absorption or emission features. 
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The spectrum of 1E1207.4-5209, on the 
contrary, appeared from the start as 
dominated by two broad absorption fea- 
tures seen, at 0.7 and 1.4 keV, both by 
Chandra (Sanwal et al, 2002) and XMM- 
Newton (Mereghetti et al, 2002). To bet- 
ter understand the nature of such fea- 
tures, XMM-Newton devoted two orbits to 
1E1207. 4-5209, for a total observing time 
of 257,303 sec. In the two MOS EPIC cam- 
eras the source yielded 74,600 and 76,700 
photons in the energy range 0.2 - 3.5 keV, 
while the pn camera recorded 208,000 pho- 
tons, time-tagged to allow for timing stud- 
ies. Analysis of this long observation, while 
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confirming the two phase-dependent ab- 
sorption lines at 0.7 and 1.4 keV, unveiled 
a statistically significant third line at ^2.1 
keV, as well as a possible fourth feature at 
2.8 keV. The nearly perfect 1:2:3:4 ratio 
of the line centroids, as well as the phase 
variation, naturally following the pulsar B- 
field rotation, strongly suggest that such 
lines are due to cyclotron resonance scat- 
tering (Bignami et al. 2003). A year after 
the XMM-Newton observation, IE 1207.4- 
5209 was deeply scrutinized by Chandra for 
^300 ksec, thus logging an overall effective 
exposure slightly larger than that devoted 
to the most popular NSs so far, such as 
Crab, Vela and RX J1856.5-3754. 

The release of a new, and improved, in- 
strument calibration software prompted us 
to revisit the XMM-Newton data set. Here 
we shall briefly report on some highlights of 
our new spectral and temporal analysis (see 
De Luca et al. 2004 for details). In addition, 
the "best" X-ray positional information is 
used, in conjunction with a new optical ob- 
servation, to investigate the source optical 
behaviour. 

2. Timing Analysis 

After converting the arrival times of the 
208,000 pn photons to the Solar System 
Barycenter, we searched the period range 
from 424.12 to 424.14 ms using both a fold- 
ing algorithm with 8 phase bins and the 
Rayleigh test. The best period value and 
its uncertainty (P = 424.13076±0.00002 
ms) were determined following the proce- 
dure outlined in Mereghetti et al. (2002). 
Comparing the new period measurement 
of IE 1207.4-5209 with that obtained 
with Chandra in January 2000 (Pavlov et 
al. 2002b), we obtain a period derivative 
P=(1.4±0.3)xl0~ 14 s s _1 . However, Fig.Q] 
(left panel) shows that the P value rests 
totally on the first Chandra period mea- 
surement. Using only the 3 most recent val- 
ues, the period derivative is unconstrained. 
Thus, we cannot exclude that the observed 
spin-down, based on only a few sparse 
measurements, be affected by glitches or 



Doppler shifts induced by orbital motion. 
Questioning the object's P would have far 
reaching consequences for the understand- 
ing of 1E1207.4-5209 since the serious dis- 
crepancy between the pulsar characteritic 
age (t c ~ 5 x 10 5 yrs) and the SNR age 
(tsnr ~ 7 kyrs) is entirely based on the 
value inferred from the measurements sum- 
marized in Fig. ^ (left panel) . 

To study the energy dependence of the 
pulse profile, we divided the data in four 
channels with approximately 52,000 counts 
each: 0.2-0.52 keV, 0.52-0.82 keV, 0.82- 
1.14 keV and 1.14-3.5 keV. The pulse pro- 
files in the different energy ranges (Fig. ^ 
right panel) show a broad, nearly sinusoidal 
shape, with a pulsed fraction varying from 
~ 3 to ~ 11 % in the four energy intervals. 
It is worth noting that the minimum pulsed 
fraction is found in the 0.20-0.52 keV en- 
ergy range, the only portion of the spec- 
trum free from absorption lines. Indeed, 
Fig- ID (right panel) is an independent con- 
firmation of the findings of Bignami et al 
(2003) who ascribed the source pulsation 
to the absorption lines phase variation. 

Finally, comparing the shapes of the 
light curves of Fig. ^ (right panel), we see 
for the first time a phase shift of nearly 
90° between the profile in the lowest en- 
ergy range (<0.52 keV) and those at higher 
energies. 

3. Spectral analysis 

The improved understanding of the in- 
struments (concerning, e.g., its Quantum 
Efficiency, its Charge Transfer Inefficiency, 
and its redistribution function) imple- 
mented in the most recent SAS release used 
here (SAS v5.4.1) yields rather significant 
differences in the low energy (E<0.6 keV) 
part of the spectrum with respect to our 
previous analysis (Bignami et al. 2003). 
This allowed for an update of the best fit- 
ting parameters for both the continuum 
and the lines. The reader is referred to De 
Luca et al. (2004) for a complete discus- 
sion of the spectral analysis of the three 
EPIC cameras. Here we shall focus on the 
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Fig. 1. (Left) Period history of IE 1207.4-5209 . Circles are the Chandra measurements, 
crosses XMM-Newton. (Right) Folded light curve of IE 1207.4-5209 in four energy ranges. 
Fitting the pulse profiles with a function consisting of a constant plus a sinusoid, we obtain 
pulsed fraction values of (2.7±0.7)% in the 0.2-0.52 keV energy range, (11.1±0.7)% for 
0.52-0.82 keV , (7.1±0.7)% for 0.82-1.14 keV, and (7.1±0.7)% 1.14-3.5 keV. 



pn data set, which contributes more than 
half of the source statistics. The best fit- 
ting continuum model comes from the sum 
of two blackbody functions. The cooler one 
has a temperature kT=0.163 ± 0.003 keV 
and an emitting radius R=4.6 ±0.1 km, 
while the hotter one has kT=0.319 ± 0.002 
and R=0.83±0.03 km. Four absorption fea- 
tures are clearly seen in the spectrum of 
IE 1207.4-5209 (see FigEJ) at the harmon- 
ically spaced energies of ^0.7 keV, ~1.4 
keV, ~2.1 keV and -2.8 keV. The dif- 
ferent spectral continuum model, resulting 
from the improved calibrations, yields also 
a more significant detection of the third and 
fourth features with respect to our previ- 
ous analysis. Using a simple gaussian in ab- 
sorption, we estimated with an F-test that 
the 2.1 keV and the 2.8 keV features have 
a chance occurrence probability of ~10 -9 
and ~ 10~ 3 , respectively. 

As a further step, we have studied the 
variations of the spectrum of IE 1207.4- 



5209 as a function of the pulse phase. 
Following Bignami et al. (2003), we se- 
lected the phase intervals corresponding to 
the peak (phase interval 0.40-0.65 with re- 
spect to Fig. ^ right panel), the declin- 
ing part (phase 0.65-0.90), the minimum 
(phase 0.90-1.15) and the rising part (phase 
0.15-0.40) of the folded light curve. The re- 
sulting spectra (see Fig. OJ) were fitted al- 
lowing both the thermal continuum and the 
lines to vary. The main results can be sum- 
marized as follows: 

— The two components of the continuum 
vary slightly both in temperature and in 
flux with the pulse phase, accounting at 
most for ^3-5% of the source pulsation. 
An inspection of Fig. 03 clearly shows 
that the phase variation of the features 
is largely responsible for the observed 
pulsation of the source. 

— the features are strongly dependent on 
the pulse phase and show significant 
variations in width, depth and shape 
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Fig. 2. The spectrum of IE 1207.4-5209 as seen by the pn camera. The upper panel 
shows the data points together with the best fitting continuum model folded with the 
instrumental response. The lower panel shows the residuals in units of standard deviations 
from the best fitting continuum. The presence of four absorption features at ~ 0.7, ~ 1.4, 
-2.1,-2.8 keV is evident. 



The central energy of the 0.7 keV fea- 
ture varies by 6% at most, while dis- 
placements of the other features are not 
significant. The 2.8 keV feature is only 
marginally detected during the mini- 
mum and the rise intervals. 
— The relative intensity of the first three 
features varies with the phase. 

4. Optimizing the X-ray position 

To derive the sky coordinates of IE 1207.4- 
5209, we computed independently for the 
MOS1 and MOS2 cameras the boresight 
correction to be applied to the default 
EPIC astrometry. We used the Guide Star 
Catalog II (GSC-II 1 ) to select, amongst our 
—200 serendipitous detections, 6 sources 
with a stellar counterpart to be used to 
correct the EPIC astrometry. The rms er- 

1 http: / /www- 
gsss.stsci.edu/gsc/gsc2/GSC2home.htm 



ror between the refined X-ray and GSC- 
II positions is —1 arcsec per coordinate. 
The resulting MOS1 position of IE 1207.4- 
5209 is a j 2 ooo = 12/il0m00.91s, (5j2ooo = 
— 52°26'28.8" with an overall error radius of 
1.5 arcsec. The MOS2 position is aj2ooo = 
12/il0ra00.84s, <5 J20 oo = -52°26'27.6" , also 
with an uncertainty of 1.5 arcsec, fully con- 
sistent with the MOS1 coordinates. 

In order to obtain an independent mea- 
surement on the position of IE 1207.4-5209, 
we have retrieved from the Chandra Data 
Archive a public dataset relative to a recent 
(June 2003) ACIS observation (20 ksec in 
imaging mode) of the target. As in the case 
of the EPIC data, we used the positions 
of the serendipitous sources in the field to 
refine the astrometry (see De Luca et al. 
2004 for details). The best Chandra/ ACIS 
position of IE 1207.4-5209 is aj 2 ooo = 
12W0m00.826s, <5 J20 oo = -52°26'28.43" 
with an uncertainty of 0.6". The more accu- 
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Fig. 3. (Left) Comparison of the four phase-resolved spectra extracted from pn data. 
Colours indicate the phase-interval: black, minimum; red, rise; green, peak; blue, decline. 
The peak of the total light curve corresponds to the phase interval where the absorp- 
tion features are at their minimum, while the light curve minimum happens when the 
absorption features are more important. (Right) Residuals in units of sigma obtained by 
comparing the data with the best fit continuum model for the phase-resolved spectra. 
The pulse phase variations in width, depth and shape of the absorption features are 
evident. 



rate Chandra position lies inside the inter- 
section of the MOS1 and MOS2 error cir- 
cles, confirming the correctness of the anal- 
ysis and the absence of systematics. 

5. Search for the optical counterpart 

The field of IE 1207-5209 was observed 
with the 8.2-mcter UT-1 Telescope (Antu) 
of the ESO VLT (Paranal Observatory). 
Observations were performed with the 
the FOcal Reducer and Spectrograph 1 
(FORS1) instrument. Images were acquired 
through the Bessel V and R filters for a 
total integration time of ~2 and 3 hrs, 
respectively. Fig. 0] shows the inner por- 
tion of the combined FORS1 1^-band im- 
age centered on the target position, with 
the ACIS, MOS1 and MOS2 error circles 
superimposed. A faint object (marked with 
the two ticks in Fig.fJ} is detected just out- 



side the southern edge of the MOS1 error 
circle. It showed variability during the time 
span covered by our observations. Also, its 
position falls more than 2 arcsec away from 
the most probable region. For both reasons, 
we can rule it out as a potential counterpart 
of IE 1207.4-5209 . 

No candidate counterpart is detected 
within the Chandra error circle (nor in 
the intersection of the MOS ones) down 
to R~27.1 and V~27.3, which we assume 
as upper limits on the optical flux of IE 
1207—5209. For an X-ray derived interstel- 
lar absorption of Ay = 0.65 at b~ 10° and 
a distance of 2 kpc, our VLT u.l. rules out 
any hypothetical "normal" stellar compan- 
ion (M>15 both in R and in V). If we as- 
sume that IE 1207-5209 is indeed isolated, 
we can derive a neutron star optical lumi- 
nosity < 3.4 1 28 erg s" 1 or < 4.6 10~ 6 of 
its rotational energy loss, a value similar to 
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Fig. 4. F0RS1 F-band image centered on 
the target position. The error circles from 
Chandra/ACIS (0.6" radius), MOS1 and 
MOS2 (2 arcsec radius) have been super- 
imposed. 



those of middle-aged INSs. Since the VLT 
flux upper limits are > 100 higher than the 
extrapolation of the XMM-Newton black- 
body (De Luca et al. 2004), they do not 
constrain our object's optical spectrum. 

6. Conclusions 

IE 1207.4-5209 is a radio-quiet NS at the 
center of a well known SNR. However, 
the longer we observe such an object, the 
harder it becomes to make it fit into a 
known class of the NS family. The long ob- 
servations devoted to IE 1207.4-5209 both 
by Chandra and by XMM-Newton have 
unveiled a number of unique and some- 
what contradictory characteristics that, at 
the moment, defy standard theoretical in- 
terpretations. 

The first and outmost problem seems 
to be the interpretation of the period evo- 
lution. The usual, simplistic, monotonic fit 
yields a P=(1.4±0.3) x 10~ 14 s s" 1 that, 
combined with the period value, makes 



IE 1207.4-5209 much too old (> 50 times) 
for its SNR. Moreover, the value of the 
star's magnetic field inferred on the ba- 
sis of the classical dipole braking formula 
(B=2.6±0.3xl0 12 G) is significantly differ- 
ent from the value obtained from the cy- 
clotron absorption lines interpreted both in 
terms of electrons (B— 8 x 10 10 G) as well 
as protons (B- 1.6 x 10 14 G). 

More creative interpretations of the pe- 
riod evolution (e.g. Dc Luca ct al. 2004, 
Zavlin et al. 2004) require scenarios such 
as bynary system or peculiar glitching. 
However, such types of behaviour are al- 
ready tightly constrained by the lack of an 
optical counterpart. 

Are the unique characteristics of 
IE 1207.4-5209 pointing towards a new 
class of NSs? 

References 

Bignami, G.F., Caraveo, P.A., De Luca, A., 

Mereghetti, S. 2003 Nature, 423, 72 
Becker W. and Aschenbach, B. 2002, in 

Neutron Stars, Pulsars, and Supernova 

Remnants, W. Becker et al. eds., p. 64 
De Luca, A., Mereghetti, S., Caraveo, P.A., 

Moroni, M., Mignani, R.P., Bignami, 

G.F. 2004, A&A, in press 
Mereghetti, S., De Luca, A., Caraveo, P.A., 

Becker, W., Mignani, R., Bignami, G.F. 

2002, ApJ, 581, 1280 
Pavlov, G.G., Zavlin, V.E., Sanwal, D. 

2002a, in Neutron Stars, Pulsars, and 

Supernova Remnants, W. Becker et al. 

eds., p.273 

Pavlov G.G., Zavlin, V.E., Sanwal, D., 
Truemper, J. 2002b, ApJL, 569, L95 

Sanwal D., Pavlov G.G., Zavlin V.E & 
Teter M.A. 2002, ApJL, 574, L61 

Zavlin V.E, Pavlov G.G., Sanwal D. 2004, 
submitted to ApJ, astro-ph/0312096 



